Mon. Not. R. Astron. Soc. 000, [I]-?? (2002) Printed 5 February 2008 (MN IATrX style file v2.2) 



A Census of the Carina Nebula. I: Cumulative Energy 
Input from Massive Stars 



Nathan Smith* 

Center for Astrophysics and Space Astronomy, University of Colorado, 389 UCB, Boulder, CO 80309, USA 



Accepted 0000, Received 0000, in original form 0000 



ABSTRACT 

The Carina Nebula (NGC 3372) is our richest nearby laboratory in which to study 
feedback through UV radiation and stellar winds from very massive stars during the 
formation of an OB association, at an early phase in the evolution of the surround- 
ing proto-superbubble before supernova explosions have influenced the environment. 
This feedback is triggering successive generations of new star formation around the 
periphery of the nebula, while simultaneously evaporating the gas and dust reservoirs 
out of which young stars are trying to accrete material. This paper takes inventory of 
the combined effect from all the known massive stars that power the Carina Nebula 
through their total ionizing flux and integrated mechanical energy from their stellar 
winds. Carina is close enough and accessible enough that spectral types for individual 
stars are available, and many close binary and multiple systems have recently been 
spatially resolved, so that one can simply add them. Adopting values from the litera- 
ture for corresponding spectral types, the present-day total ionizing photon luminosity 
produced by the 65 O stars and 3 WNL stars in Carina is Qh ~10 51 s _1 , the total 
bolometric luminosity of all stars earlier than B2 is 2.5 xlO 7 Lq, and the total me- 
chanical luminosity of stellar winds is L$w — 10 5 Lq. The total Qh was about 25% 
higher when rj Carinae was on the main sequence, before it and its companion were 
surrounded by its obscuring dust shell; for the first 3 Myr, the net ionizing flux of the 
70 O stars in Carina was about 150 times greater than in the Orion Nebula. About 
400-500 Mq has been contributed to the H n region by stellar wind mass loss during 
the past 3 Myr. Values for Qh and Lsw are also given for the individual clusters 
Trl4, 15, and 16, and BolO and 11, which are more relevant on smaller spatial scales 
than the total values for the whole nebula. 

Key words: H II regions — ISM: bubbles — ISM: individual (NGC 3372) — stars: 
early-type — stars: formation — stars: mass-loss — stars: winds, outflows 



1 INTRODUCTION 

Feedback from young massive stars may play an integral 
role in star and planet formation. Most stars are born in 
the vicinity of the hottest and most massive stars spawned 
only from giant molecular clouds, and the effects of feed- 
back from these massive stars cannot be studied in nearby 
quiescent regions of star formation. Even the nearest H II 
region, the Orion Nebula, may not be representative of the 
extreme environments where most stars are born, since it is 
dominated by just a single 06 dwarf. 

A much more extreme collection of stars can be found in 
the Carina Nebula (NGC 3372). Its star clusters are almost 
as spectacular as those of 30 Dor in the LMC (Massey & 
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Hunter 1998), or objects in our own Galaxy like the Arches 
cluster near the Galactic center (Najarro et al. 2004; Figer et 
al. 1999, 2002) and NGC 3603 (Moffat et al. 2002). However, 
these other regions are too distant for detailed studies of 
small-scale phenomena like irradiated protoplanetary disks 
and jets, and their study is hampered by considerably more 
extinction. The low extinction toward Carina combined with 
its proximity and rich nebular content provide a worthwhile 
trade-off for its less concentrated star clusters. While this 
may be a liability for investigating the upper end of the 
initial mass function, the fact that Carina's stars are more of 
a loose aggregate may be an advantage for studying feedback 
and triggered star formation. Carina provides a snapshot 
of an OB association in the making, which may be more 
representative of the environments in which most stars form 
than are bound super star clusters. It may also provide an 



2 N. Smith 



early analog of the local Galactic environment that evolved 
into Gould's Belt. 

Recent studies suggest that at the present epoch there 
exists an upper mass limit to the most massive stars of 
~150 M (Figer 2005, Weidner & Kroupa 2004, Kroupa 
& Weidner 2005; Oey & Clarke 2005). While the star clus- 
ters that power Carina may not be as dense as the Arches 
cluster or NGC 3603, Carina does contain several examples 
of stars that probably began their lives with 100-150 Mq. 
Among these are rj Carinae, the prototypical 02 supergiant 
HD 93129A, three late-type hydrogen-rich Wolf-Rayet stars 
(WNL stars), and the remaining original members of the 
03 spectral class for which this spectral type was first in- 
troduced (Walborn 1973; although see Walborn et al. 2002). 
For these colossal stars, lifetimes of ~3 Myr before they ex- 
plode as supernovae (SNe) are shorter than the time it takes 
to clear away their natal molecular material. Consequently, 
we already have a situation in Carina where the most mas- 
sive members like rj Car are approaching their imminent 
demise while new stars are being born from dense molecu- 
lar gas only 5-20 pc away. In the next 1-2 Myr, there will 
be several very energetic SNe in the Carina Nebula, which 
will carve out an even larger cavity in the ISM and form a 
giant superbubble in the Galactic plane, and may pollute 
protoplanetary disks with nuclear-processed ejecta. In the 
mean time, studying this rich region provides a snapshot 
of a young proto-superbubble (Smith et al. 2000) energized 
only by UV radiation and stellar winds, immediately before 
the disruptive SN-dominated phase. 

In just the past decade, the Carina Nebula has been 
recognized as a hotbed of ongoing active star formation. It 
provides a laboratory to study several star formation phe- 
nomena in great detail, all of which have recently been iden- 
tified in this region: (1) evaporating protoplanetary disks 
(the so-called "proplyds"), small cometary clouds, or glob- 
ules (Smith et al. 2003a; Smith 2002&; Brooks et al. 2000, 
2006; Cox & Bronfman 1995), (2) the erosion of large dust 
pillars and the triggering of a second generation of embed- 
ded star formation within them (Smith et al. 2000, 2005; 
Rathborne et al. 2004; Megeath et al. 1996), (3) irradiated 
Herbig-Haro (HH) jets that are a signpost of active accretion 
(Smith et al. 2004), (4) photodissociation regions (PDRs) on 
the surfaces of molecular clouds across the region (Brooks 
et al. 2003, 2006, 1998; Rathborne et al. 2002; Smith et al. 
2000; Mizutani et al. 2004), and on the largest scales, (5) 
the early formation of a fledgeling superbubble (Smith et al. 
2000). These phenomena trace a second generation of stars, 
perhaps triggered by feedback from the first generation clus- 
ters. Carina is near enough that an upcoming program with 
the Hubble Space Telescope will undertake the same types of 
detailed studies of feedback that have been done in Orion 
at visual and near-IR wavelengths, but sampling a more ex- 
treme environment. 

All of these phenomena respond to the cumulative UV 
radiation and/or stellar winds from star clusters that power 
the H II region, and a census of this energy input is essential 
in order to understand the relationship between them. Es- 
timating the total stellar energy input in Carina is a more 
difficult problem than in smaller H II regions like Orion, 
where a single star dominates the photoevaporation of pro- 
plyds and the irradiation of HH jets. Instead, Carina has 
dozens of massive O stars over many parsecs. Although the 



remarkable stellar content of the Carina Nebula has been 
discussed many times in the literature (e.g., Walborn 1995, 
2002, 2005; Feinstein 1995; Massey & Johnson 1993), there is 
no complete and up-to-date census of the massive stars, and 
especially there has been no estimate of their collective UV 
radiation and mechanical (wind) luminosity that powers the 
region. This is the main purpose of the present investigation. 
A second paper in this series will compare the total energy 
input from stars with the apparent energy budget inferred 
from observations of the surrounding Carina Nebula. 



2 ADOPTED PARAMETERS 

The basic strategy in this paper will be to use the observed 
spectral type of each massive star (spectral types as late as 
B2 will be considered, although their influence is minimal) 
to assign a corresponding value for the hydrogen-ionizing 
photon luminosity Qh, as well as the mass- loss rate M and 
terminal wind speed Voo, and hence, the wind mechanical 
luminosity Lsw=(l/2)MVl . Table 1 lists the values that 
will be adopted for stars later than 03. In Table 1, values of 
log L and Qh for each spectral type are taken from the "ob- 
served" temperature-scale models of Martins, Schaerer, & 
Hillier (2005). Values for M are taken from Repolust, Puis, 
& Herrero (2004), and terminal velocities for each spectral 
type are taken by a rough fit to the average values listed by 
Prinja, Barlow, & Howarth (1990). 1 

The stellar content of the Carina Nebula is truly re- 
markable and accessible, and several of the more spectac- 
ular stars have received enough attention that the desired 
parameters have been derived individually for them. In some 
cases, when noted below, individual studies may supersede 
the average values adopted in Table 1. Many of these stars 
are, in fact, benchmarks for their spectral types in the three 
papers just mentioned. In addition, the three WNL stars 
in Carina are of particular interest, as are HD 93129A and 
t] Car (when it was on the main sequence), since these are 
the most luminous sources that have the most influence. The 
specific stellar content of the nebula is described in more de- 
tail in the following section. 



3 STELLAR CONTENT 

Tables 2-6 compile the stellar content powering the Carina 
Nebula. References for the stellar content of each cluster are 
noted in each table. Several reviews exist in the literature 
(especially Walborn 1995, 2005; Feinstein 1995); these were 
invaluable in compiling the tables. Every effort has been 
made to provide the most recent information with regard 
to binarity /multiplicity through spectral types or through 
high spatial resolution techniques, especially for Trl4 and 
16. This, of course, does not preclude the possibility that 
additional companions may be identified, and their influence 
should be considered by readers in that case. 

Following Walborn (1995) and Tapia et al. (1988), it 

Given the spread in terminal velocities listed by Prinja et al. 
1990, we adopted the same velocities for all three luminosity 
classes in Table 1. This is a reasonable approximation for the 
earliest spectral types that dominate the energy budget. 
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Table 1. Adopted Parameters for Each Spectral Type 



Spect. 


Lum. 


log L 


log Qh 


M 




L SW 


Type 


Class 


(L ) 


(s L ) 


CIO" 6 M yr- 1 ) 


/. —Is 

(km s 1 ) 


(L ) 


03 


V 


5.84 


49.64 


3.1 


3160 


2450 


03.5 




5.76 


49.54 


2.5 


3080 


1870 


04 


V 


5.67 


49.44 


2.0 


2990 


1410 


04.5 


V 


5.58 


49.33 


1.6 


2870 


1040 


OS 


V 


5.49 


49.22 


1.3 


2760 


780 


OS. 5 


V 


5.41 


49.10 


1.0 


2650 


550 


06 


V 


5.32 


48.99 


0.8 


2560 


410 


06.5 


V 


5.23 


48.88 


0.6 


2460 


290 


07 


V 


5.14 


48.75 


0.4 


2320 


170 


07. 5 


V 


5.05 


48.61 


0.35 


2210 


135 


OS 


V 


4.96 


48.44 


0.26 


2100 


91 


08.5 


V 


4.86 


48.27 


0.22 


1960 


67 


09 


V 


4.77 


48.06 


0.18 


1820 


47 


09. 5 


V 


4.68 


47.88 


0.15 


1660 


33 


BO 


V 


4.57 


47.70 


0.12 


1510 


22 


BO. 5 


V 


4.47 


47.50 


0.10 


1460 


17 


Bl 


V 


4.37 


47.28 


0.08 


1180 


8.8 


B1.5 


V 


4.28 


47.05 


0.06 


960 


4.4 


B2 


V 


4.19 


46.80 


0.05 


750 


2.2 


03 


III 


5.96 


49.77 


6.5 


3160 


5130 


03.5 




5.91 


49.71 


5.5 


3080 


4120 


C)4 


III 


5.85 




4.7 


2990 


3320 


04.5 


III 


5.79 


49.56 


4.0 


2870 


2600 


OS 


III 


5.73 


49.48 


3.4 


2760 


2050 


OS. 5 


III 


5.67 


49.40 


2.7 


2650 


1500 


06 


III 


5.61 


49.32 


2.2 


2560 


1140 


06.5 


III 


5.54 


49.23 


1.8 


2460 


860 


07 


III 


5.48 


49.13 


1.5 


2320 


640 


07. 5 


III 


5.42 


49.01 


1.2 


2210 


460 


OS 


III 


5.35 


48.88 


0.8 


2100 


280 


08.5 


III 


5.28 


48.75 


0.7 


1960 


210 


09 


III 


5.21 


48.65 


0.5 


1820 


130 


09. 5 


III 


5.15 


48.42 


0.4 


1660 


87 


BO 


III 


5.08 


48.28 


0.3 


1510 


54 


BO. 5 


III 


5.00 


48.10 


0.25 


1460 


42 


Bl 


III 


4.93 


47.90 


0.20 


1180 


22 


B1.5 


III 


4.86 


47.68 


0.15 


960 


11 


B2 


III 


4.78 


47.44 


0.11 


750 


4.9 


03 


I 


5.99 


49.78 


13.0 


3160 


10300 


03. 5 


I 


5.96 


49.74 


11.3 


3080 


8470 


04 


I 


5.93 


49.70 


10.0 


2990 


7060 


04.5 


I 


5.90 


49.66 


8.9 


2870 


5790 


OS 


I 


5.87 


49.62 


8.0 


2760 


4810 


05. 5 


I 


5.84 


49.58 


7.1 


2650 


3940 


06 


I 


5.81 


49.52 


6.4 


2560 


3310 


06.5 


I 


5.78 


49.46 


5.7 


2460 


2730 


07 


I 


5.75 


49.41 


5.0 


2320 


2130 


07. 5 


I 


7.72 


49.31 


4.3 


2210 


1660 


08 


I 


5.68 


49.25 


3.6 


1960 


1090 


08.5 


I 


5.65 


49.19 


3.0 


1820 


790 


09 


I 


5.61 


49.11 


2.5 


1660 


540 


09. 5 


I 


5.57 


49.00 


1.9 


1510 


340 



References. — L, Q H : Martins ct al. (2005); Voo : Prinja ct al. (1990); M: 
Repolust et al. (2004). L and Qh for B stars were extrapolated from Martins 
et al. (2005) following Crowthcr (2005). 



is assumed that all these clusters share a common distance 
to the Carina Nebula, taken here to be 2.3±0.1 kpc (Smith 
2002a; Walborn 1995; Allen & Hillier 1993). The distance of 
2.3 kpc determined for n Carinae itself is very reliable, hav- 
ing been derived from proper motions and Doppler velocities 
of the expanding Homunculus (Smith 2002a; Allen & Hillier 
1993). Furthermore, r\ Car is known to be at the same dis- 
tance as the Keyhole nebula, since its reflected spectrum is 
seen across the Keyhole (Walborn & Liller 1977; Lopez & 
Meaburn 1986). Differences in photometric distance moduli 
given in the literature are probably due to anomalous or 
varying reddening laws, patchy extinction, small age differ- 
ences, or other factors. For example, photometric distances 
as high as 3.7 kpc have been given for Trl5 (Walborn 1973), 



requiring that Trl5 be several hundred pc behind, the rest of 
the Carina nebula. This seems implausible, although not im- 
possible, for several reasons, including the fact that Trl5 is 
seen along the same line of sight as dense gas and dust struc- 
tures, some of which are dust pillars that point toward Trl6. 
Yet, stars in Trl5 have an average E(B — V)~0.5 (Feinstein 
ct al. 1980) - roughly the same as Trl4 and 16 (Walborn 
1995; Feinstein et al. 1973) - which could not be the case if 
Trl5 were far behind the rest of the nebula and obscured. 
Similarly, relatively large distances to Trl4 (Vazquez et al. 
1996) are implausible. Even though Trl4 must be somewhat 
more distant than Trl6 (perhaps as much as 20 pc) because 
silhouette objects projected in front of it point toward Trl6 
(Smith et al. 2003a), it cannot be several hundred parsecs 
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Content o 


f Trumpler 


16 and C 


ollinc 


Icr 228 






Name 


Sp. 


log L 


log Qh 




M 




Lsw 








V b 1 


do~ 6 

V-IA' 


o y > 


(km s~ 1 ) 


y^Q) 


rj Car-now 


LBV 


6.67 






1000 


600 


28700 


■q Car-MS 


02 If (?) 


6.57 


50.36 




8.32 


3200 


6780 


r] Car-B 


OS V (?) 


5.49 


49.2 




1.3 


2760 


780 


HD 93162 


WN6ha 


6.22 


50.0 




10.5 


2480 


5140 


HD 93131 


WN6ha 


5.94 


49.7 




13.8 


2160 


5130 


HD 92740A 


WN7ha 


6.01 


49.8 




14.5 


1785 


3680 


HD 92740B 


08-9 V 


4.86 


48.27 




0.22 


1960 


67 


HD 93205A 


03.5 V((f)) 


5.76 


49.54 




2.5 


3080 


1870 


HD 93205B 


08 V 


4.96 


48.44 




0.26 


2100 


91 


HD 93250 


03.5 V((f)) 


5.76 


49.54 




2.5 


3080 


1870 


HDE 303308 


04 V((f)) 


5.67 


49.44 




2.0 


2990 


1410 


HD 93204 


OS V((f)) 


5.49 


49.22 




1.3 


2760 


780 


Trl6-244 


04 If 


5.93 


49.70 




10.0 


2990 


7060 


CPD -59 2600 


06 V((f)) 


5.32 


48.99 




0.8 


2560 


410 


CPD -59 2603A 


07 V((f)) 


5.14 


48.75 




0.4 


2320 


170 


CPD -59 2603B 


09.5 V 


4.68 


47.88 




0.15 


1660 


33 


CPD -59 2603C 


B0. 2 V 


4.57 


47.70 




0.12 


1510 


22 


CPD -59 2628A 


09.5 V 


4.68 


47.88 




0.15 


1660 


33 


CPD -59 2628B 


B0. 3 V 


4.57 


47.70 




0.12 


1510 


22 


CPD -59 2635A 


08 V 


4.96 


48.44 




0.26 


2100 


91 


CPD -59 2635B 


09.5 V 


4.68 


47.88 




0.15 


1660 


33 


CPD -59 2636A 


07 V 


5.14 


48.75 




0.4 


2320 


170 


CPD -59 2636B 


08 V 


4.96 


48.44 




0.26 


2100 


91 


CPD -59 2636C 


09 V 


4.77 


48.06 




0.18 


1820 


47 


CPD -59 2641 


OS V 


5.49 


49.22 




1.3 


2760 


780 


HDE 303311 


OS V 


5.49 


49.22 




1.3 


2760 


780 


HDE 305536 


09 V 


4.77 


48.06 




0.18 


1820 


47 


HD 93027 


09.5 V 


4.68 


47.88 




0.15 


1660 


33 


HD 93028 


09 V 


4.77 


48.06 




0.18 


1820 


47 


HD 93130 


06 111(f) 


5.61 


49.32 




2.2 


2560 


1140 


HD 93146 


06.5 V((f)) 


5.23 


48.88 




0.6 


2460 


290 


HD 93160 


06 III((f)) 


5.61 


49.32 




2.2 


2560 


1140 


HD 93161Aa 


OS V 


4.96 


48.44 




0.26 


2100 


91 


HD 93161Ab 


09 V 


4.77 


48.06 




0.18 


1820 


47 


HD 93161B 


06.5 V((f)) 


5.23 


48.88 




0.6 


2460 


290 


HD 93206A 


09.7 Ib:(n) 


5.57 


49.00 




1.9 


1510 


340 


HD 93206B 


09 III 


5.21 


48.65 




0.5 


1820 


130 


HD 93222 


07 111(f) 


5.48 


49.13 




1.5 


2320 


640 


HD 93343 


07 V(n) 


5.14 


48.75 




0.4 


2320 


170 


HD 93403A 


OS III(f)vai 


5.73 


49.48 




3.4 


2760 


2050 


HD 93403B 


07 V 


5.14 


48.75 




0.4 


2320 


170 


Trl6-2 


Bl V 


4.37 


47.28 




0.08 


1180 


8.8 


Trl6-3 


08.5 V 


4.86 


48.27 




0.22 


1960 


67 


Trl6-4 


Bl V 


4.37 


47.28 




0.08 


1180 


8.8 


Trl6-5 


Bl V 


4.37 


47.28 




0.08 


1180 


8.8 


Trl6-8 


B0. 5 V 


4.47 


47.50 




0.10 


1460 


17 


Trl6-9 


09.5 V 


4.68 


47.88 




0.15 


1660 


33 


Trl6-10 


B0 V 


4.57 


47.70 




0.12 


1510 


22 


Trl6-ll 


B1.5 V 


4.28 


47.05 




0.06 


960 


4.4 


Trl6-12 


Bl V 


4.37 


47.28 




0.08 


1180 


8.8 


Trl6-13 


Bl V 


4.37 


47.28 




0.08 


1180 


8.8 


Trl6-14 


B0. 5 V 


4.47 


47.50 




0.10 


1460 


17 


Trl6-16 


Bl V 


4.37 


47.28 




0.08 


1180 


8.8 


Trl6-17 


Bl V 


4.37 


47.28 




0.08 


1180 


8.8 


Trl6-19 


09.5 V 


4.68 


47.88 




0.15 


1660 


33 


Trl6-20 


Bl V 


4.37 


47.28 




0.08 


1180 


8.8 


Trl6-21 


08 V 


4.96 


48.44 




0.26 


2100 


91 


Trl6-22 


08.5 V 


4.86 


48.27 




0.22 


1960 


67 


Trl6-23 


07 V 


5.14 


48.75 




0.4 


2320 


170 


Trl6-24 


B2 V 


4.19 


46.80 




0.05 


750 


2.2 


Trl6-27 


Bl V ? 


4.37 


47.28 




0.08 


1180 


8.8 


Trl6-28 


B2 V 


4.19 


46.80 




0.05 


750 


2.2 


Trl6-29 


B2 V 


4.19 


46.80 




0.05 


750 


2.2 


Trl6-31 


B0. 5 V 


4.47 


47.50 




0.10 


1460 


17 


Trl6-37 


Bl V ? 


4.37 


47.28 




0.08 


1180 


8.8 


Trl6-74 


Bl V 


4.37 


47.28 




0.08 


1180 


8.8 


Trl6-76 


B2 V 


4.19 


46.80 




0.05 


750 


2.2 


Trl6-94 


B1.5 V 


4.28 


47.05 




0.06 


960 


4.4 


Trl6-11B 


08.5 V 


4.86 


48.27 




0.22 


1960 


67 


Trl6-122 


B1.5 V 


4.28 


47.05 




0.06 


960 


4.4 


Trl6-124 


Bl V 


4.37 


47.28 




0.08 


1180 


8.8 


Trl6-126 


09 V 


4.77 


48.06 




0.18 


1820 


47 


Trl6-127 


09 V 


4.77 


48.06 




0.18 


1820 


47 


Trl6-245 


B0 V 


4.57 


47.70 




0.12 


1510 


22 



Note — rj Car's hypothetical value for Qh on the main sequence is taken by adopting the 
same Q H /^Bol ra ti° as f° r HD 93129Aa. References. — Nclan ct al. (2004); Masscy Sz Johnson 
(1993); Walborn (1973, 1995; 2002, 2005); Penny ct al. (1993); Repolust ct al. (2004); Naze ct 
al. (2005; for HD93161); Smith ct al. (1996); Crowthcr ct al. (1995, 2002; although Q H and M 
values for the WNL stars arc from Crowthcr 2005, private comm., modified to include clumping 
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farther away. The outermost filaments of the Carina region 
span a distance of ~120 pc (Smith et al. 2000), yet Trl4 and 
16 are projected only about 5 pc from each other on the sky, 
and even the Keyhole Nebula near Trl6 shows signs of ion- 
ization from the northwest (i.e. from Trl4). 

Additional clusters that are sometimes associated with 
Carina, such as NGC 3324 and NGC 3293, are ignored in 
the present analysis. They do not appear to lie within the 
nebular structures that define the outermost boundaries of 
the Carina Nebula, and the stars of NGC 3324 are associated 
with their own separate circular H II region. 



3.1 Trumpler 16 

Energy input in Carina is dominated by the massive cluster 
Tr 16, of which the infamous massive star rj Car is the most 
luminous member. This is clear even before a detailed look 
at its stellar content, since optical images of Carina show 
numerous sculpted dust pillars and other elongated struc- 
tures located across several degrees of the sky, and nearly 
all of them point back toward rj Car and the center of Trl6 
(Smith et al. 2000, 2003a, 2004, 2005; Walborn 2002). 

As noted by Walborn (1995) and others, Cr 228 and Cr 
232 are probably part of Trl6, where their apparent separa- 
tion on the sky is an artifact caused by dust lanes on the near 
side of the nebula. Thus, Table 2 includes Cr 228 and 232 
as members of a larger Trl6 cluster. Table 2 also includes 
the 05 V star HDE 303311, which is part of Cr232, as a 
member of Trl6, even though its membership is somewhat 
uncertain (Walborn 1995). In this and other cases, whether 
this star is considered part of Trl6 or Trl4 has no net effect 
on the global energetics of the region, but it may inflate the 
relative importance of Trl6. 

Also included within Trl6 in Table 2 are the three WNL 
stars in the region. HD 93162 is almost certainly a member 
of Trl6. The two other WNL stars, HD 92740 and HD 93131, 
are quite a bit farther away from rj Car on the sky. Walborn 
(1995) notes that they are both roughly equidistant (~20 
pc) from rj Car, and that they could easily have traveled 
this far in 2 Myr with a kick of ~10 km s _1 . Previously, on 
the main-sequence, they presumably had 02 spectral types 
like HD 93129A, and somewhat different Qh values. In order 
to account for this stellar evolution, the cumulative "main- 
sequence" values in Table 7 (see §4) were calculated using 
Qh values lowered by 0.1 dex for each WNL star. This 
is roughly correct if they were ~0.1 dex less luminous on 
the main-sequence (Maeder & Meynet 2000) , with the same 
Qh/LboI ratio as HD93129A. 

The age of Trl6 is probably about 2-3 Myr. This comes 
from the fact that it contains H-rich WNL stars but no He- 
rich WR stars, and that rj Car has a likely ZAMS mass >100 
Mq and has already evolved off the main-sequence — as well 
as the fact that there is no clear evidence for a supernova 
having exploded yet. 

When considering the cumulative effect of the massive 
stars in Trl6, we must bear in mind the peculiar present 
state of its most luminous member, rj Car. As we observe it 
today, rj Car is a pitiful source of UV photons for the larger 
region. Whatever ionizing radiation is able to escape its very 



dense optically-thick stellar wind 2 (Hillier et al. 2001; Smith 
et al. 2003&) is eventually absorbed by dust in the surround- 
ing Homunculus nebula and converted into infrared radia- 
tion (Smith et al. 2003c; Cox et al. 1995). However, more 
than 10 4 years ago, before it entered its presently-observed 
luminous blue variable (LBV) phase and was on the main 
sequence — or even more than 160 years ago before the Ho- 
munculus existed — rj Car was a much stronger source of UV 
radiation. Indeed, it probably contributed ~20% of the to- 
tal ionizing photons for the entire region (see below) . By the 
same token, its main-sequence stellar wind was weaker than 
its presently-observed LBV wind. Its spectral type prior to 
the LBV phase is of course unknown, but given its extreme 
luminosity, it is reasonable to assume that it was an 02 su- 
pergiant like HD 93129A while on the main sequence. From 
stellar evolution tracks for very massive stars (e.g., Young 
2005; Maeder & Meynet 2000), it is likely that rj Car was 
~0.1 dex less luminous at that time. The value for Qh at 
that time is then calculated assuming this lower Lso!, and 
assuming the same ratio of Qh/LboI as for HD93129A. 

It is also likely that rj Car has a close companion star in 
a 5.5 yr eccentric orbit (Damineli et al. 2000; Corcoran 2005; 
Duncan et al. 1995). While its UV radiation does not escape 
the Homunculus at the present time, it did contribute to the 
ionization of the larger region throughout the 2-3 Myr life- 
time of Trl6. The nature of rj Car's companion star is highly 
uncertain and a matter of current debate (Pittard & Corco- 
ran 2002; Ishibashi et al. 1999; Davidson 1999; Duncan et al. 
1995); a tentative spectral type of 05 V is adopted in Table 
2. The bolometric luminosity of rj Car (the primary star) in 
Table 2 has the contribution of its companion subtracted. 

3.2 Trumpler 14 

Trumpler 14 is a smaller, much more compact, and per- 
haps somewhat younger cluster than Trl6. Some members 
of Trl6, like HD 93161, 93160, and 93250, are projected 
nearby on the sky (Walborn 1973) and are often mistaken 
as members of Trl4. As noted above, it is likely that Trl4 
is perhaps 10-20 pc more distant than Trl6, since silhou- 
ette objects seen in front of it point toward Trl6 (Smith et 
al. 2003), but it is still part of the central engine that pow- 
ers the Carina Nebula. Its most conspicuous member is the 
prototypical 02 If* supergiant HD 93129A. 

Because of its compact size, its relatively faint upper 
main sequence, and certain spectral diagnostics, Trl4 may 
be 1-2 Myr younger than Trl6 (Walborn 1995). Penny et 
al. (1993) estimate ages younger than 1 Myr for HD 93129A 
and HD 93128, while Vazquez et al. (1996) find a cluster age 
of ~1.5 Myr. In addition, Trl4 remains closer to its natal 
molecular cloud than Trl6 (e.g., Brooks et al. 2003). 

The stellar content summarized in Table 3 is compiled 
from several references listed at the end of the table, with 
particular attention to recent results on the binarity and 
multiplicity of stars in the cluster (see Nelan et al. 2004; 
Walborn 2005), including a new companion to HD 93129A 
(Nolan et al. 2004) . The greatest source of uncertainty in the 
net output of this cluster is in the parameters for HD93129A 

2 The mass-loss rate for rj Car listed in Table 2 is from Hillier et 
al. (2001), and includes the effects of clumping in the wind. 
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Table 3. Stellar Content of Trumpler 14 



Name 


Sp. 


log L 


log Qh 


M 








Type 


( L ©) 


(s ) 


(10~ 6 M yr- 1 ) 


(km s 1 ) 


(L Q ) 


HD 93129Aa 


02 If* 


6.17 


49.96 


8.32 


3200 


6780 


HD 93129Ab 


03.5 V 


5.76 


49.54 


2.5 


3080 


1870 


HD 93129B 


03.5 V((f+)) 


5.76 


49.54 


2.5 


3080 


1870 


HD 93128 


03.5 V((f+)) 


5.76 


49.54 


2.5 


3080 


1870 


CPD-58 2611 


06 V((f)) 


5.32 


48.99 


0.8 


2560 


410 


CPD-58 2620 


06.5 V((f)) 


5.23 


48.88 


0.6 


2460 


290 


Trl4-3 


BO. 5 IV- V 


4.50 


47.70 


0.15 


1460 


26 


Trl4-4 


BO V 


4.57 


47.70 


0.12 


1510 


22 


Trl4-5 


09 V 


4.77 


48.06 


0.18 


1820 


47 


Trl4-6 


Bl V 


4.37 


47.28 


0.08 


1180 


8.8 


Trl4-9 


08 V 


4.96 


48.44 


0.26 


2100 


91 


Trl4-18 


BO V 


4.57 


47.70 


0.12 


1510 


22 


Trl4-21 


09 V 


4.77 


48.06 


0.18 


1820 


47 


Trl4-27 


09 V 


4.77 


48.06 


0.18 


1820 


47 


Trl4-30 


BO III-IV 


4.90 


48.10 


0.25 


1510 


45 



References. — Nelan et al. (2004); Massey & Johnson (1993); Walborn (1973, 1995; 2005); Penny et al. 
(1993); Repolust et al. (2004). 



Table 4. Stellar Content of Trumpler 15 



Name 


Sp. 


log L 


log Qh 


M 


Voo 






Type 


(L©) 




(10- 6 Mq yr" 1 ) 


(km s- 1 ) 


(L©) 


HDE 303304 


08 V 


4.96 


48.44 


0.26 


2100 


91 


HD 93249 


08 II/09 III 


5.28 


48.75 


0.7 


1960 


210 


HD 93342 


09 III 


5.21 


48.65 


0.5 


1820 


130 


Trl5-18 


09 I-II 


5.61 


49.11 


2.5 


1660 


540 


Trl5-2 


09 III 


5.21 


48.65 


0.5 


1820 


130 


Trl5-19 


09 V 


4.77 


48.06 


0.18 


1820 


47 


HD 93190 


B0 IV 


4.83 


48.00 


0.21 


1510 


38 


Trl5-3 


B2 V 


4.19 


46.80 


0.05 


750 


2.2 


Trl5-4 


Bl V 


4.37 


47.28 


0.08 


1180 


8.8 


Trl5-7 


B2 V 


4.19 


46.80 


0.05 


750 


2.2 


Trl5-9 


Bl V 


4.37 


47.28 


0.08 


1180 


8.8 


Trl5-10 


B2 V 


4.19 


46.80 


0.05 


750 


2.2 


Trl5-13 


Bl V 


4.37 


47.28 


0.08 


1180 


8.8 


Trl5-14 


B2 V 


4.19 


46.80 


0.05 


750 


2.2 


Trl5-15 


B0. 5 IV-V 


4.47 


47.50 


0.10 


1460 


17 


Trl5-20 


Bl V ? 


4.37 


47.28 


0.08 


1180 


8.8 


Trl5-21 


B0 III 


5.08 


48.28 


0.3 


1510 


54 


Trl5-23 


B0 V 


4.57 


47.70 


0.12 


1510 


22 


Trl5-26 


B0. 5 V 


4.47 


47.50 


0.10 


1460 


17 


References. — 


Fcinstcin et al. 


(1980); Morrcll et al. 


(1988). 







itself. The value of log Ls !=6.17 is adopted from Repolust 
ct al. (2004), and the value for Qh was calculated from 
the luminosity using the revised T e // calibration scale of 
Martins et al. (2005). 



3.3 Trumpler 15 

Trl5 is an open cluster located about 20' north of Trl6 on 
the sky. It is thought to be at the same distance and has 
roughly the same reddening as Trl6, but it is somewhat 
older, with a likely age of 6±3 Myr (Feinstein et al. 1980; 
Morrell et al. 1988). It is seen near a very bright red star 
on the sky - the M2 la supergiant RT Car, which is proba- 
bly not a member of the cluster. Table 4 adopts the cluster 
members listed by Feinstein et al. (1980) and Morrell et al. 
(1988), including the 08 V star HD 303304 which was listed 
as a possible nonmember by Feinstein et al. (1980), as well 
the 09 giant HD 93342 and the highly-reddened 09 I-II star 
Trl5-18, which were listed as possible members. Table 4 does 
not include the non-member 05 III star HD 93403, which is 
taken to be a member of Trl6. 



3.4 Bochum 10 

BolO is a relatively meager cluster located roughly 40' north- 
west of Trl6, with a handful of OB stars and a probable age 
of ~7 Myr (Feinstein 1981; Fitzgerald et al. 1987). Thus, it 
may be significantly older than Tr 14 and 16, but perhaps 
not much older than Trl5. Feinstein found a large distance 
of 3.6 pc, but a later analysis by Fitzgerald et al. (1987) fa- 
vored a closer distance of 2.5 kpc, placing it inside the Carina 
complex associated with some nearby nebulosity. Because of 
its larger age of ~7 Myr, it is likely that the nearby WR star 
HD 92809 may be associated with BolO, or at least with the 
larger Carina complex. Therefore, HD 92809 is included in 
Table 5 along with the other members of BolO. For the total 
FUV luminosity of HD92809 in Table 7 (see §4), we adopt 
L FUV /L Bo i=0.39 (for 912-3650 A) following the study of 
Smartt et al. (2001). 



3.5 Bochum 11 

Boll is a loose open cluster in the southeastern Carina Neb- 
ula. Fitzgerald et al. (1987) find an age for Boll of less than 
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Table 5. Stellar Content of Bochum 10 



Name 


bp. 


log L 


log Qh 


M 


Voo 


L SW 




Type 


(L©) 


(s- 1 ) 


(10- 6 M Q yr- 1 ) 


(km s" 1 ) 


(L©) 


HD 92809 


WC6 


5.32 


49.20 


16 


2280 


6620 


HD 92725 


09.5 III 


5.15 


48.42 


0.4 


1660 


87 


HD 92964 


B2 la 


5.30 


48.35 


0.5 


750 


22 


HD 92759 


BO III 


5.08 


48.28 


0.3 


1510 


54 


HD 92894 


BO IV 


4.83 


48.00 


0.21 


1510 


38 


HD 93002 


B0 IV 


4.83 


48.00 


0.21 


1510 


38 


HDE 302989 


B2 V 


4.19 


46.80 


0.05 


750 


2.2 


HDE 303190 


Bl V 


4.37 


47.28 


0.08 


1180 


8.8 


HDE 303296 


Bl V 


4.37 


47.28 


0.08 


1180 


8.8 


HDE 303297 


Bl V 


4.37 


47.28 


0.08 


1180 


8.8 


BolO-2 


B0. 5 IV 


4.74 


47.80 


0.18 


1460 


30 


BolO-7 


B2 IV 


4.49 


47.12 


0.08 


750 


3.5 


BolO-15 


Bl V 


4.37 


47.28 


0.08 


1180 


8.8 



References. — Fcinstein (1981); Fitzgerald et al. (1987); Prinja et al. (1990); Smartt et al. (2001). 

Table 6. Stellar Content of Bochum 11 



Name Sp. log L log Q H M Voo Lsw 

Type (L ) (s- 1 ) (1()- 6 M yr" 1 ) (kms- 1 ) (L ) 



HD 93632 


04-5 111(f) 


5.79 


49.56 


4.0 


2870 


2600 


Boll-2 


09 V 


4.77 


48.06 


0.18 


1820 


47 


Boll-3 


B2 V 


4.19 


46.80 


0.05 


750 


2.2 


Boll-4 


B0. 5 V 


4.47 


47.50 


0.10 


1460 


17 


Boll-5 


09 V 


4.77 


48.06 


0.18 


1820 


47 


Boll-9 


09 IV 


4.99 


48.36 


0.34 


1820 


89 


Boll-10 


09 IV 


4.99 


48.36 


0.34 


1820 


89 



References. — Fitzgerald et al. (1987); Walborn (1973). 



3 Myr. They also find that the cluster includes several pre- 
main-sequence stars, consistent with an age as young as 0.3 
Myr. This very young age is interesting in light of the fact 
that Boll is located amid what is currently the most ac- 
tive region of star formation in Carina - the so-called South 
Pillars (Smith et al. 2000; Rathborne et al. 2004). The age 
and stellar content is similar to the nearby embedded clus- 
ter known as the Treasure Chest (Smith et al. 2005; Hagele 
et al. 2004). 3 There is considerable discrepancy as to the 
spectral type of the most luminous star in Bo 11, which is 
HD 93632 or Boll-1. Fitzgerald et al. (1987) list this star as 
08 lap, whereas Walborn (1973) classifies it as 04-5 111(f). 
Walborn (1973) noted that the spectral type was somewhat 
variable, while Fitzgerald et al. (1987) noted a few prob- 
lems with their photometric classification of stars in Boll. 
A spectral type of 04-5 III(f ) is adopted in Table 6. 



4 CUMULATIVE EFFECTS 

Table 7 lists the total luminosity, ionizing flux, total FUV 
(Balmer continuum) luminosity, mass loss, and mechanical 
luminosity for each cluster, as well as the cumulative total 
of all clusters for each parameter. For the FUV luminosity 
(in all cases except where noted above), Lftjv ~0.5Lb o i is 
assumed, which is adequate considering the inherent uncer- 
tainty in the bolometric luminosity as a function of spectral 
type. Brooks et al. (2003) found a similar result for the FUV 
luminosity of Trl4 alone. More detailed stellar stmosphere 

3 Although it is not discussed in detail here, the 09.5 V star 
CPD-59°2661 in the Treasure Chest cluster is included in Table 
7. 



models typically have Lfuv /^Boi values of 0.39-0.85 for 
912-3650 A (e.g., Smartt et al. 2001; Crowther et al. 2002; 
Crowther 2005, private comm.). Three cases are given for 
Trl6 and for the total of all clusters: 

The first case correponds to the history of Carina up un- 
til recent times, when rj Car and its companion were not sur- 
rounded by a dust shell and did not have a dense LBV wind 
choking off the Lyman continuum flux, and when the WNL 
stars in Trl6 were 02 stars. For this first case, there was a 
total of 70 O-type stars in Carina, producing Qh=L15x 10 51 
s _1 . This would be the appropriate number to adopt when 
considering the history and formation of the nebula, the life- 
times of evaporating proplyds, globules, and dust pillars, 
triggered star formation by radiative driven implosion, and 
the growth of the cavity that will blow out of the Galactic 
plane as a bipolar superbubble. 

The second case is the LBV phase of i) Car during the 
last ~10 4 yr, when r) Car had a dense wind that reprocessed 
most of the Lyman continuum radiation into the Balmer 
continuum (FUV). During this phase, r\ Car was not sur- 
rounded by a dense dust shell, allowing its FUV luminosity 
to escape. Since the UV radiation of its companion could 
also escape, and since the WNL stars were no longer on 
the main sequence, the total number of O stars is listed as 
66 in Table 7. This brief phase may be important in un- 
derstanding the recent photoevaporation of neutral globules 
and photodissociation regions throughout the nebula, since 
7] Car was such an incredibly luminous FUV source. 

The third case represents the currently-observed state 
of the Carina Nebula, when rj Car and its companion are 
surrounded by an obscuring dust shell, blocking all their 
contribution to the total ionizing flux and FUV luminos- 
ity. With 7] Car's companion blocked by the Homunculus, 
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Table 7. Total Stellar Energy Input 



Cluster 


Number of 


log L 


log Qh 


log L(FUV) 


M 






O stars 


(L©) 


(s- 1 ) 


(L©) 


(lO- 6 Mq yr- 1 ) 


(L©) 


Trl6 (MS) 


47 


7.215 


50.91 


6.91 


91 


45400 


Trl6 (LBV) 


43 


7.240 


50.78 


7.05 


1083 


67000 


Trl6 (now) 


42 


7.240 


50.77 


6.79 


1083 


67000 


Trl4 


10 


6.61 


50.34 


6.31 


18.7 


13500 


Trl5 


6 


6.18 


49.56 


5.88 


5.9 


1300 


BolO 


1 


6.00 


49.42 


5.69 


18.3 


7120 


Boll 


5 


6.00 


49.64 


5.70 


5.2 


2900 


CPD-59°2661 


1 


4.68 


47.88 


4.38 


0.15 


33 


TOTAL (MS) 


70 


7.38 


51.06 


7.08 


139 


70200 


TOTAL (LBV) 


66 


7.40 


50.97 


7.18 


1131 


91800 


TOTAL (now) 


65 


7.40 


50.96 


7.00 


1131 


91800 



the effective number of O stars is reduced to 65. For this 
third case, the cumulative ionization source is Qh=9xW 50 
s _1 . This would be the appropriate number to adopt when 
considering the current UV flux incident upon evaporating 
proplyds, globules, and irradiated jets over most of the neb- 
ula. 

For most of its lifetime, gas and dust in the Carina 
Nebula has been exposed to an ionizing luminosity about 150 
times stronger than that of the Orion nebula. This changed 
160 years ago when r) Car ejected a thick dust shell that 
cut off essentially all of its UV output, and the total Qh of 
the Carina Nebula dropped by about 20% due to the loss of 
ionization from the region's most luminous member. 4 

This change illustrates the dominant role of the most 
luminous star in any region, and highlights a truly unique 
property of the Carina Nebula — it offers a laboratory where 
we can address variability in Qh over short timescales. This 
behavior has not been witnessed anywhere else in our Galaxy 
to a similar degree. We can observe some of the effects of 
this variable UV output — for example, neutral cometary 
globules and dust pillars that point at r\ Car are seen only 
in silhouette today, but their shape suggests that they were 
formed by a strong UV flux from rj Car in the past (Smith 
et al. 2003a). It is interesting to note that the outer edges 
of the Carina Nebula, beyond a radius of 160 light yr or 50 
pc (1.5xl0 20 cm) from r/ Car, have not yet seen the Great 
Eruption of 77 Car, and still see its pre-eruption UV out- 
put. Also, if the complex nested circumstellar ejecta around 
77 Car are a fair indication, this type of UV cutoff has prob- 
ably happened multiple times in the past due to previous 
outbursts, and may yet happen again. 

The main uncertainty in deriving the pre-eruption UV 
output of Carina is in the spectral type and ionizing luminos- 
ity of r\ Car itself, which has a likely uncertainty of perhaps 
±15% (roughly 3% of the total for the region) if the models 
at lower luminosity and later spectral types are correct. In 
terms of the cumulative effect of the entire stellar popula- 
tion, then, dwarf stars later than 09 can be safely ignored, 

4 Actually, the change might not have been so sudden, since 
much of the Lyman continuum flux from i) Car may have already 
dropped when it developed a thick stellar wind upon entering the 
LBV evolutionary phase in the last 10 4 yr. During this time, the 
majority of r\ Car's radiation escaped as FUV photons. 




1 10 50 

Radius (pc) 

Figure 1. Ionizing photon flux $ as a function of radius. The two 
solid curves are for Case 1 and 3 (i.e. 77 Car on the main sequence 
and at the present time). 



since their collective contribution is less than the uncertainty 
in the value for Qh from the most luminous members. For 
the other luminous members, in Case 1 (when 77 Car was on 
the main sequence), about 8% of the total Qh came from 
the 02 supergiant HD 93129Aa alone, 15% came from the 3 
WNL stars, and 15% from the remaining five 03.5 V stars. 
For case 3, corresponding to the present state of affairs, these 
contribution have changed to about 10% of the total from 
HD 93129Aa alone, 23% from the 3 WNL stars, and 19% 
from the remaining five 03.5V stars. 

Figure 1 shows the ionizing photon flux &h as a func- 
tion of radius. The two solid curves are for Case 1 and 3 (i.e. 
77 Car on the main sequence and at the present time). Be- 
yond a radius of roughly 10 pc from the center of the nebula, 
the ionizing source can be treated roughly as a point source, 
so that the photon flux drops as r~ 2 . Within a few parsecs 
of any cluster, however, the local radiation from that clus- 
ter becomes more important than the diminished collective 
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1 10 50 

Radius (pc) 

Figure 2. Same as figure 1, but for FUV radiation, plotted in 
units of the Habing flux, Go=1.6xlO — 3 ergs s _1 cm~ 2 . The three 
cases from Table 7 arc shown. 



effects from all clusters. This is due to the fact that the stel- 
lar population in Carina is a loose aggregate, rather than a 
single dense cluster. The specific case of Trl4 is shown in 
Figure 1, but within a few parsecs of any individual cluster, 
its own ionization will dominate, while the influence of more 
distant clusters diminishes. 

Figure 2 shows the radial dependence of the FUV ra- 
diation field in the Carina Nebula, compared to the Habing 
flux, G =1.6xl0" 3 ergs s" 1 cm" 2 (Habing 1968). The three 
curves correspond to r\ Car on the main sequence, in the 
LBV phase, and in its present state obscured by its own dust 
shell. The strongest FUV fluxes in the region arise when the 
LBV phase of r\ Car is considered, when most of its bolomet- 
ric luminosity escaped in the Balmer continuum. Through- 
out the Carina Nebula, this resulted in a local FUV flux that 
was roughly 75% higher than the present state when r\ Car 
is surrounded by dust. Thus, for extremely massive stars, 
their changing UV output in various evolutionary states can 
profoundly effect the radiation field of their environments, 
and consequently, the evaporation of protoplanetary disks 
and the erosion of dust pillars. Again, however, individual 
clusters dominate over the collective radiation field within a 
few parsecs of any cluster core. 

The present total mechanical luminosity from stellar 
winds is just under 10 5 Lq, or 7xl0 4 Lq during the main- 
sequence phase, having contributed a total energy of about 
2.6 x 10 52 ergs so far during the 3 Myr lifetime of the region. 
This is higher than a previous estimate of several x 10 51 ergs 
(Walborn 1982) because it includes the influence of many 
more stars. The estimated kinetic energy of the expanding 
proto-superbubble (Smith et al. 2000; Walborn 1982) is only 
5-30% of this mechanical energy from stellar winds, suggest- 
ing that most of the energy injected by stellar winds is ra- 
diated away, even in the early phases of a bubble's growth. 
The total energy radiated by early-type stars during the 



same time was over 350 times more than the mechanical en- 
ergy. This situation will change dramatically when the most 
massive stars begin to explode as supernovae in the next few 
10 5 -10 6 yr, each injecting an additional 10 51 -10 52 ergs in a 
very short time. 

With a total mass-loss rate of M « 1.4 x 10~ 4 M 
yr -1 , the mass ejected by stellar winds over the past 3 Myr is 
about 420 Mq. This is only about 10% of the mass of ionized 
gas estimated to be filling the interior of the H II region (e.g., 
Walborn 2005). Thus, the majority of ionized gas inside the 
H II region cavity probably comes from photoevaporative 
flows (e.g., Bertoldi 1989) off the surfaces of molecular clouds 
and evaporating globules embedded within the region. Note 
that the present value of M for r] Car's brief LBV wind 
phase is much greater than the cumulative mass loss of rest 
of the stars in Carina combined. 



5 COMPARISON WITH OTHER MASSIVE 
CLUSTERS 

Finally, it is interesting to see how the stellar energy input 
in Carina stacks up to other young massive star clusters. 
As noted earlier, in general, Carina's overall stellar con- 
tent is similar to the most well-studied massive clusters like 
NGC3603, Arches, and 30 Dor, except that these clusters 
are more compact than Carina. For example, Crowther & 
Dessart (1998) list log Q w =51.19 and L sw =1.5xl0 5 L Q 
for NGC3603. The total energy input in Carina is essen- 
tially the same as these values, because Crowther & Dessart 
used an older T e // scale, leading to higher Qh values for O 
stars, and did not account for wind clumping in their anal- 
ysis; if those same assumption were made here for Carina, 
the results would be identical to NGC3603 within the un- 
certainties. Carina has about 25-30% of the total ionizing 
flux of the more massive Arches cluster in the Galactic cen- 
ter (Figer et al. 2002), and roughly 25% of both the total 
Qh and Lsw of R136 in 30 Doradus given by Crowther & 
Dessart (1998). Given the similar cumulative energy input 
and stellar content, it is interesting that Carina's clusters are 
less tightly bound than the others, especially if this gives rise 
to different effects on the nearby star- forming environment. 
A future paper will investigate how the cumulative stellar 
energy input assessed directly from the observed massive 
stars in Carina compares with more indirect tracers of the 
energy budget of the larger Carina Nebula. 
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